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Abstract 

The sensitivity of calorimetric energy-energy correlation function to the medium- 
induced energy loss of fast partons in high multiplicity heavy ion interactions is demon- 
strated at the appropriate selection of events for the analysis, namely, the availability 
of one high-pT jet in an event at least and using the procedure of "thermal" background 
subtraction. Without jet trigger this correlation function manifests the global structure of 
transverse energy flux: the correlator is isotropic for central collisions, and for non-central 
collisions it is sensitive to the azimuthal anisotropy of energy flow reproducing its Fourier 
harmonics but with the coefficients squared. 
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1 Introduction 



The QGP creation and properties of this new type of superdense matter are extensively dis- 
cussed in the current literature [1] dedicated to high energy heavy ion physics. In the last several 
years many new phenomena have been observed at RHIC, such as strong elliptic flow [2] and 
suppression of the high transverse momentum two particle back-to-back correlations [3]. These 
observations together with other important evidences support the idea that a dense partonic 
matter (QGP) has been created in such high energy nuclear collisions. 

The center of mass energy for heavy ion collisions at the LHC will exceed that at RHIC 
by a factor of about 30. This provides exciting opportunities for addressing unique physics 
issues in a new energy domain, where hard and semi-hard QCD multi-particle production can 
certainly dominate over underlying soft events [4]. The methodological advantage of azimuthal 
jet observables is that one needs to reconstruct only the azimuthal position of the jet in heavy ion 
event, that can be done with fine enough resolution (slightly worse as compared with pp case but 
still lesser than the typical azimuthal size of a calorimeter tower [4, 5]). However the observation 
of the jet azimuthal anisotropy due to medium-induced partonic energy loss [6, 7] requires event- 
by-event determination of the nuclear event plane. Recently the ability to reconstruct the 
reaction plane using calorimetric measurements has been shown in [8] as well as the possibility 
to observe jet azimuthal anisotropy without direct determination of the event plane, considering 
the second [9] and higher [10] order correlators between the azimuthal position of the jet axis 
and the angles of the particles not incorporated in the jet. This calorimetric information allows 
us to study also the energy-energy correlations which are sensitive to partonic energy loss (as 
it has been pointed out in [11]) and this investigation does not demand also the reconstruction 
of the event plane which is still far from to be completed. 

The observed at RHIC suppression of the high transverse momentum two particle back- 
to-back correlations [3] supports strongly the believe that the calorimetric measurements of 
energy-energy correlations at LHC will be able to provide a new additional information about 
medium-induced energy loss of fast partons. Namely these correlations are the subject of the 
present investigation. 

2 Energy-energy correlation functions in e + e~, hadronic 
and nuclear collisions 

The energy-energy correlation function S has been used by all for LEP experiments [12] at 
CERN and the SLD experiment [13] at SLAC to measure the strong coupling constant a s in 
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e + e -annihilation at the Z° resonance with a high accuracy. £ is defined as a function of the 
angle x between two particles % and j in the following form 

rfcos(x) Acos( X )iV cvcntc ^ t ^ E* ' 1 J 

where E is the total energy of the event, and £j are the energies of the particles i and j. 
The sum runs over all pairs i, j with cos(x) in a bin width Acos(x): 

cos(x) - Acos(x)/2 < cos(x) < cos(x) + Acos(x)/2. 

a is the total cross section for e + e~ — > hadrons. The limits A cos(x) — > and N cvcnt — > oo have 
to be taken in (1). 

This function can be calculated in perturbative QCD as a series in a s : 

1 dE(x) a s{lA ,/ x , { a s (^)\ 2 ( a , / A* 



cr d cos(x) 27T 



^ (A,ln(|)A( X ) + S( X ))+0(a2), (2) 



where /? = (33 — 2n/)/3, is the number of active flavours at the energy E 1 . The first order 
term has been calculated by Basham et al. [14] from the well-known one gluon emission 
diagrams 7*, Z° — > gg^ with the result 

A( x ) = C F (1 + cu) 3 ^^((2 - 6u 2 ) ln(l + 1/w) + 6cu - 3), (3) 

where Cp = 4/3, = cot 2 (x)/2 and x is the angle between any of the three partons. This 
allowed one to determine the strong coupling constant directly from a fit to the energy-energy 
correlations in e + e~-annihilation, since £ is proportional to a s in the first order. 

In hadronic and nuclear collisions jets are produced by hard scattering of partons. In this 
case it is convenient to introduce the transverse energy-energy correlations which depend very 
weakly on the structure functions [15] and manifest directly the topology of events. Thus, for 
instance, in high-p T two-jet events the correlation function is peaked at the azimuthal angle 
ip — 0° and tp = 180°, while for the isotropic background this function is independent of ip. 

Utilization of hard jet characteristics to investigate QGP in heavy ion collisions is com- 
plicated because of a huge multiplicity of "thermal" secondary particles in an event. Various 
estimations give from 1500 to 6000 charged particles per rapidity unit in a central Pb-Pb col- 
lision at LHC energy, and jets can be really reconstructed against the background of energy 
flux beginning from some threshold jet transverse energy E^ ~ 50 — 100 GeV [4, 5, 16]. The 
transverse energy-energy correlation function is just sensitive to the jet quenching as well as 
to the global structure of energy flux (i.e. its anisotropy for non-central collisions) if we select 
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events by an appropriate way. The generalization of (1) for calorimetric measurements of the 
energy flow is straightforward: 

dT, T (<p) _ 1 E Ti E T ( i+n ) 

d<p ~ A<p N event (E™y ' [ ) 

where 

i 

is the total transverse energy in N calorimetric sectors covering the full azimuth, E?i is the 
transverse energy deposition in a sector % (i = 1, N) in the considered pseudo-rapidity region 
1 77 1, n = [tp/Aip] (the integer part of the number ip/A<p), A<p ~ 0.1 is the typical azimuthal size 
of a calorimetric sector. In the continuous limit Ap — > 0, N = [2n/A(p\ — > 00 the equation (4) 
reads 

dE T ((p) 1 ^ 1 2 f dEf cIEt , , \ , cN 



dp iV event e ^ t (£f ls ) 2 / # ^ ^ d<P 

/dEx dEj- 
d<p — — ((f)), — — is the distribution of transverse energy flow over azimuthal 
dcp dcp 

angle 0. 

The events with high-p-r jet production are rare events and their contribution to £t is 
negligible without the special jet trigger. As a result, the transverse energy-energy correlator is 
expected to be independent of ip for central collisions due to the isotropy of energy flow on the 
whole. For non-central collisions with the clearly visible elliptic flow of the transverse energy, 

^ = ^[1 + 2^^(2(0-^))], (6) 

this correlator dT,r/dip is independent of the reaction plane angle i/jr and is calculated explicitly: 

^> = £p+ *«-<*)]. (7) 

The strength of collective flow (6) being determined by the value of v 2 , and the intensity of 
oscillations in (7) - - by v\ ! Moreover one can prove that dY> T /d(p reproduces all Fourier 
harmonics of transverse energy flow decomposition but with the coefficients squared. 

3 The model to simulate heavy ion events without and 
with jets at the LHC 

In order to see the manifestation of the jet structure we must cut the background events, i.e. con- 
sider the sum in (4) over only the events containing at least one jet with > i£ T (threshold) ~ 
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100 GeV and subtract the background from soft thermal particles using its isotropy on the 
whole. We demonstrate the productivity and effectiveness of such strategy in the framework 
of our well worked-out model of jet passing through a medium applied early to the calculation 
of various observables sensitive to the partonic energy loss: the impact parameter dependence 
of jet production [17], the mono/dijet rate enhancement [18, 19], the dijet rate dependence 
on the angular jet cone [20], the elliptic coefficient of jet azimuthal anisotropy [6, 7, 8, 9, 10], 
the anti-correlation between softening jet fragmentation function and suppression of the jet 
rate [21]. For model details one can refer to these mentioned papers (mainly, [17, 20, 21]). 
Here we note only the main steps essential for the present investigation. 

PYTHIA_6.2 [22] was used to generate the initial jet distributions in nucleon-nucleon sub- 
collisions at y/s = 5.5 TeV. After that event-by-event Monte-Carlo simulation of rescattering 
and energy loss of jet partons in QGP was performed. The approach relies on accumulative 
energy losses, when gluon radiation is associated with each scattering in expanding medium 
together including the interference effect by the modified radiation spectrum as a function of 
decreasing temperature. Such numerical simulation of free path of a hard jet in QGP allows 
any kinematical characteristic distributions of jets in the final state to be obtained. Besides the 
different scenarios of medium evolution can be considered. In each i-th scattering a fast parton 
loses energy collisionally and radiatively, Aej = tj/(2m ) + iq^ where the transfer momentum 
squared tj is simulated according to the differential cross section for elastic scattering of a 
parton with energy E off the "thermal" partons with energy (or effective mass) m ~ 3T <C E 
at temperature T, and Ui is simulated according to the energy spectrum of coherent medium- 
induced gluon radiation in BDMS formalism [23]. Finally we suppose that in every event the 
energy of an initial parton decreases by the value AE = J2i Ae;. 

The medium was treated as a boost-invariant longitudinally expanding quark-gluon fluid, 
and partons as being produced on a hyper-surface of equal proper times r [24]. For certainty 
we used the initial conditions for the gluon-dominated plasma formation expected for central 
Pb— Pb collisions at LHC [25]: t ~ 0.1 fm/c, T ~ 1 GeV. For non-central collisions we 
suggest the proportionality of the initial energy density to the ratio of nuclear overlap function 
and effective transverse area of nuclear overlapping [17]. 

The energy-energy correlator depends on not only an absolute value of partonic energy 
loss, but also on an angular spectrum of in-medium radiated gluons. Since coherent Landau- 
Pomeranchuk-Migdal radiation induces a strong dependence of the radiative energy loss of a jet 
on the angular cone size [20, 23, 26, 27, 28], it will soften particle energy distributions inside the 
jet, increase the multiplicity of secondary particles, and to a lesser degree, affect the total jet 
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energy. On the other hand, collisional energy loss turns out to be practically independent of jet 
cone size and causes the loss of total jet energy, because the bulk of "thermal" particles knocked 
out of the dense matter by elastic scatterings fly away in almost transverse direction relative 
to the jet axis [20]. Thus although the radiative energy loss of an energetic parton dominates 
over the collisional loss by up to an order of magnitude, the relative contribution of collisional 
loss of a jet growths with increasing jet cone size due to essentially different angular structure 
of loss for two mechanisms [20]. Moreover, the total energy loss of a jet will be sensitive to 
the experimental capabilities to detect low-pr particles - products of soft gluon fragmentation: 
thresholds for giving signal in calorimeters, influence of strong magnetic field, etc. [5]. 

Since the full treatment of angular spectrum of emitted gluons is rather sophisticated and 
model-dependent [20, 23, 26, 27, 28], we considered two simple parameterizations of distribu- 
tion of in-medium radiated gluons over the emission angle 9. The "small-angular" radiation 
spectrum was parameterized in the form 



where 9 ~ 5° is the typical angle of coherent gluon radiation estimated in work [20]. The 
"broad-angular" spectrum has the form 



We believe that such simplified treatment here is enough to demonstrate the sensitivity of the 
energy-energy correlator to the medium- induced partonic energy loss. 

The following kinematical cuts on the jet transverse energy and pseudo-rapidity were ap- 
plied: > 100 GeV and | 7]j ct | < 1.5. After this the dijet event is imposed upon the Pb— Pb 
event, which was generated using the fast Monte-Carlo simulation procedure [9, 16, 29] giving a 
hadron (charged and neutral pion, kaon and proton) spectrum as a superposition of the thermal 
distribution and collective flow. To be definite, we fixed the following "freeze-out" parameters: 
the temperature Tf = 140 MeV, the collective longitudinal rapidity Y^ ax = 3 and the collective 
transverse rapidity Y™ ax = 1. We set the Poisson multiplicity distribution. For non-central 
collisions, the impact parameter dependence of the multiplicity was taken into account by a 
simple way, just suggesting that the mean multiplicity of particles is proportional to the nuclear 
overlap function. We also suggested [9] that the spatial ellipticity of the "freeze-out" region 
is directly related to the initial spatial ellipticity of the nuclear overlap zone. Such "scaling" 
allows one to avoid using additional parameters and, at the same time, results in an elliptic 
anisotropy of particle and energy flow due to the dependence of the effective transverse size 




(8) 
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of the "freeze-out" region on the azimuthal angle of a "hadronic liquid" element. Obtained in 
such a way the azimuthal distribution of particles is described well by the elliptic form (6) for 
the domain of reasonable impact parameter values. 

4 Numerical results and discussion 

At first, we became convinced that for non-central collisions with the elliptic anisotropy of 
particle and energy flow (generated in the framework of simple Monte-Carlo procedure described 
above) the energy-energy correlation function dY^r/dtp (5) is closely followed the formula (7). 

Further, only central collisions (in which the effect of jet quenching is maximum and "ther- 
mal" background is azimuthally isotropic) will be considered. In order to extract jet-like energy- 
energy correlator from the "thermal" background in high multiplicity environment, the energy 
deposition in each calorimeter sector i is recalculated event-by-event as 



where E Ti (ip) = (Eft? imal (ip) + E^(ip)) is the average energy deposition in a calorimeter 



sector and D T = y (E Ti ((p)) 2 — (E Ti (ip)) 2 is the energy dispersion in the given event, N is the 
total number of calorimeter sectors. If Eft™ becomes negative, it is set to zero. Thus only the 
energy deposition higher than the confidential interval k ■ Dt around the mean value is taken 
into account decreasing the background influence by such a way. Let us emphasize that due 
to event-by-event background fluctuations, subtracting not only the average energy deposition 
but also its dispersion with some factor k > is necessary (although this results also in some 
reduction of a signal itself). The similar procedure is applying to increase the efficiency of jet 
reconstruction in heavy ion collisions at LHC [4, 5]. 

To be specific, we consider the geometry of the CMS detector [30] at LHC. The central 
("barrel") part of the CMS calorimetric system covers the pseudo-rapidity region \rj\ < 1.5, the 
segmentation of electromagnetic and hadron calorimeters being A77 x A<fr = 0.0174 x 0.0174 and 
Ar] x A0 = 0.0872 x 0.0872 respectively [30]. In order to reproduce roughly the experimental 
conditions (not including real detector effects, but just assuming calorimeter hermeticity), we 
applied (10) to the energy deposition of the generated particles, integrated over the rapidity in 
N = 72 sectors (according to the number of sectors in the hadron calorimeter) covering the full 
azimuth. Then the energy-energy correlation function (4) is calculated for the events containing 
at least one jet with E£ > _&r (threshold) = 100 GeV in the considered kinematical region. 
The final jet energy is defined here as the total transverse energy of final particles collected 



E^itp) = E Tl &) - E Ti { V ) — k ■ D T , 



(10) 
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around the direction of a leading particle inside the cone R = ^Ar/' 2 + Aip 2 = 0.5, where rj and 
if are the pseudo-rapidity and the azimuthal angle respectively. Note that the estimated event 
rate in the CMS acceptance, ~ 10 7 jets with Et > 100 GeV in a one month LHC run with 
lead beams [4, 5], will be large enough to carefully study the energy-energy correlations over 
the whole azimuthal range. 

Before discussing numerical results, let us enlarge upon the role of background fluctua- 
tions for our analysis. The simple estimation on this can be obtained from the following 
considerations. The average (over events) energy deposition in a calorimeter sector for the 
mean total particle multiplicity (N p ) and mean particle transverse momentum (pt) is ap- 
proximately (£^ crmal ) « (p T ) ■ (N p ) /N. Then the energy dispersion in a calorimeter sec- 
tor can be roughly estimated as (7 ti iermal ~ (p T } . (N p ) /N. For realistic values N = 72, 
(Pt) = 0.55 GeV/c and (N p ) ~ 20000 (corresponds to the charged particle density per unit 
rapidity dN ± /dy{y = 0) = 5000 and Ar]=3), we get (£$ ermal ) ~ 150 GeV and ^ hcrmal ~ 10 
GeV. Namely these background parameters were fixed for our calculations. Note that a re- 
ally measurable energy flux (and correspondingly absolute value of its fluctuations) can be 
essentially less than above values, because of limited experimental capabilities to detect low-p T 
particles which giving bulk of total multiplicity. In particular, in CMS most of charged low-p^ 
particles may be cleared out of the central calorimeters by the strong magnetic field. Thus the 
analysis of quality of extracting jet-like energy-energy correlator from background fluctuations 
based on the detailed simulation of detector responses (together with the optimization of jet 
detection threshold) should be performed for each specific experiment. 

The results of numerical simulations are presented on Fig.l for the cases without and with 
partonic energy loss, two parameterizations of distribution on gluon emission angles (8) and 

(9) being used. The important issue is the procedure of event-by-event background subtraction 

(10) (with the factor k = 2 here) allows the extraction of jet-like energy-energy correlation 
function to be done even in high multiplicity heavy ion events. The efficiency of such strategy 
is closely related with the following. Due to the background isotropy, the average (over events) 
"thermal" energy deposition in the given calorimeter sector, ^E^ crma1 ^, is approximately equal 
to the average "thermal" energy deposition in a sector in the given event, E^f rmal (ip), at the 
large enough N and numbers of particles in sectors. 

The correlation function is sensitive to the medium- induced partonic energy loss and angu- 
lar spectrum of gluon radiation. Three features of medium-modified energy-energy correlation 
function can be noted: 

(1) Moderate broadening of the near-side jet region if < 0.5 (more pronounced for the "small- 
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angular" radiation), which is related with discussed in the recent literature jet shape modifica- 
tion [31]. 

(2) Significant strengthening in the wide region of azimuthal angles around tp ~ 71/ 2 (more 
pronounced for the "broad-angular" radiation) at the relatively small, but hopefully still sta- 
tistically reliable signal values. 

(3) Small additional suppression of back-to-back correlations for ip ~ n (the analog of acopla- 
narity [32]) as compared to the original (in pp) suppression due to the initial and final state 
gluon radiation. 

The observation of above modifications demands both the high enough statistics (> 10 6 
events) and the fine azimuthal resolution of jet position (< 0.1 rad) which, however, are expected 
to be attainable in calorimetric measurements at LHC [4, 5]. 

Note that for events with the trigger on the pair of energ eticjets (E^(l,2) > ^(threshold) = 
100 GeV) the calculated correlator is practically insensitive to the partonic energy loss, mainly, 
because such jets must lose small amount of energy to be visible simultaneously. 

5 Conclusions 

In summary, at the special selection of events for the analysis (at least one high-p^ jet) and the 
procedure of event-by-event background subtraction in high multiplicity heavy ion collisions, the 
transverse energy-energy correlator is sensitive to the partonic energy loss and angular spectrum 
of radiated gluons. Medium-modified energy-energy correlation function manifests significant 
strengthening in the wide interval of azimuthal angles around 7r/2, moderate broadening of the 
near-side jet region ip < 0.5 and weak additional suppression of back-to-back correlations for 
ip ~ 7r. Without jet trigger this correlation function shows the global structure of transverse 
energy flux: the correlator is isotropic for central collisions and for non-central collisions it is 
sensitive to the azimuthal anisotropy of energy flow reproducing its Fourier harmonics but with 
the coefficients squared. We believe that such energy-energy correlation analysis may be useful 
at LHC data processing. 
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Figure 1: The transverse energy-energy correlator d^/dp as a function of azimuthal angle 
tp without (solid histogram) and with medium-induced partonic energy loss for the "small- 
angular" (8) (dotted histogram) and the "broad-angular" (9) (dashed histogram) parameteri- 
zations of emitted gluon spectrum in central Pb— Pb collisions. Applied kinematical cuts are 
described in the text. 
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